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Abstract 

We present analytical and numerical properties of the coherently rotating Affleck- 
Dine condensate formed along MSSM flat directions with the radiative correc- 
tions to the mass terms included. We analyse the pressure p of the condensate, 
which is known to be negative if the potential grows slower than the field squared. 
We show that ellipticity of the orbit of the rotating field also affects the pres- 
sure. For circularly orbiting field the p = while the smallest negative value, 
corresponding to the most unstable configuration, is obtained for a coherently 
oscillating field. The AD condensate is known to fragment into non-top ological 
solitons called Q-balls. We also study equilibration of Q-ball distribution re- 
sulting from the fragmentation of the condensate analytically. We find that 
equilibration is dependent on the energy-to-charge ratio x of the condensate. 
We find the allowed range of x numerically and deduce that equilibration is 
likely to happen. 
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1 Introduction 



The Affleck-Dine (AD) mechanism [p]] provides a model for generating the observed 

baryon asymmetry of the universe in the framework of supersymmetry. In this scenario 

some squarks and/or sleptons acquire a large vacuum expectation value (VEV) along 

a flat direction (AD field) of the scalar potential of the MSSM during inflationary 

epoch. A baryon/lepton number violating operator is induced by new physics at a 

high scale together with a large C and CP violating phase, which is provided by the 

initial VEV along the flat direction. Together with the out of equilibrium conditions 

after inflation the three requirements for the generation of baryon asymmetry are 

satisfied |2j . The AD field starts oscillating^ once its effective mass exceeds the Hubble 

expansion rate H . At the same time a baryon/lepton number violating operator 

produces a torque towards the phase direction which leads to a spiral motion of the 

VEV in the complex field space. This results in a baryon/lepton asymmetry in a form 

of spatially homogeneous condensate once the comoving number density of the AD 

particles is frozen at sufficiently late times. 

The dynamics after the AD condensate formation is determined by the mass of 

the flat direction. This is essentially given by the soft SUSY breaking mass terms 

of the MSSM since the Hubble induced mass term is negligible at this stage. When 

radiative corrections are taken into account the mass term typically tends to grow 

slower than quadratically [||, [|. This induces in the condensate a negative pressure 

and makes it unstable with respect to spatial perturbations. Due to the instabilities 

the AD condensate fragments into lumps of charge, which relax into Q-balls. Q-ball 

is a field configuration which minimizes the energy with non-zero charge ||. The 

fragmentation process is highly non-linear and has to be solved numerically. AD 

condensate fragmentation has been extensively studied in [7], |^, [TIJ where it was 

noted that the original charge of the AD condensate is stored in the Q-balls. It was 

also noted that the relative amount of Q-balls and anti-Q-balls is determined by the 

dimensionless ratio of energy per charge times the scalar mass, denoted as x in this 
1 Strictly speaking AD field is rotating if the condensate has a non-zero charge. Pure oscillation 
has charge zero. But in the literature the term oscillation has become standard. 
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paper. For x > 1 there appear anti-Q-balls in addition to Q-balls. If x 3> 1 the 
distribution of Q-balls approaches a thermal one [Q, |IOj. Therefore it is important to 
know the energy of the condensate in order to study the fragmentation process. 

In this paper we calculate numerically the properties of the AD condensate: charge 
density in the comoving volume, energy-to-charge ratio x and the pressure-to-energy 
ratio w which gives the equation of state. When the charge density is known we can 
estimate the resulting charge-to-entropy ratio, which is known from nucleosynthesis 
calculations. This has been studied in |TT|, [TjJ. However, x and w have not been 
calculated previously. We find that generically \x\ > 1.1 in gravity mediated SUSY 
breaking models. For gauge mediated models the situation is more involved because 
the mass is almost a constant at the high energy scale. This effectively leads to a 
larger negative pressure which makes the energy of the condensate grow rapidly. This 
effect is also seen in the gravity mediated case but is much milder. 

These same considerations apply to any extension of the MSSM. For the extensions 
there can exist flat directions which carry charges different from baryon or lepton 
number but the same considerations apply provided that their radiative corrections 
give a mass term that again grows slower than quadratically. Q-balls, which do not 



carry baryon or lepton number, could act as self-interacting dark matter |L3, 14, [15 



This paper is organized as follows. In Sec. II we review the main features of 



AD mechanism ||11|| , calculate the pressure of a rotating condensate by extending the 
analysis of Jl6| for oscillating condensates, review the relevant properties of Q-balls 
and re-analyze the reaction rates to show that for large x a thermal equilibrium is a 
generic feature. In Sec. Ill we give the details and results of the numerical analysis. 
In Sec. IV we give the conclusions. 



2 General properties of AD mechanism 
2.1 The evolution of flat directions 

The potential for the flat direction is induced by the soft SUSY breaking mass terms, 
the A-terms of the MSSM and by the finite energy SUSY breaking due to the super- 
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gravity coupling of the flat direction to the inflaton superfield in the early universe. 
All these effects result in a potential of the form (TT| 



v W = v mW - + [ Am ^t: H ^ + h,j + i$r<«>, (i) 

where ch ~ 1, M p ~ 10 18 GeV is a high mass scale which here has been chosen as the 
reduced Planck mass, d is the dimension of the non-renormalizable operator that lifts 
the flat direction (we consider the cases d — 4, 5, 6, 7) and A, A and a are in general 
complex constants 0(1). Here we make the difference between D- and F-term inflation 
]T7| : in the D-term case the Hubble induced A-term is absent, i.e. a = 0, and in the 



F-term case a ~ 1. V m (<&) is the mass term of the flat direction which, including the 
radiative corrections, in case of gravity mediated SUSY breaking ||J is 

l$l 



V m ($) = m 



2" 

2 i + mog 



3/2 



M 2 



m (2) 



where m 3 / 2 ~ 1 TeV is the gravitino mass, K ~ —0.1 ... — 0.01 when considering un- 
stable flat directions and M ~ (M^- 3 m 3/2 /| A|) 1/(d ~ 2) is the associated renormalization 



scale [18]. In the case of gauge mediated SUSY breaking 



1$ 



2 s 



V m ($)=m%log 1 + Hr » ( 3 ) 



m 



where m$ ~ 1 . . . 100 TeV; in this case the gravitino mass m 3 / 2 < 1 GeV. One should 
note that in order to have an AD condensate formation the parameters of the A-term 
are constrained: in the gravity mediated case \A\ < d and in the gauge mediated case 
|^4| m 3/2 ^ 10 -4 , 10 _7 m$ when d — 4, 6. If the condition in the gauge mediated case 
is not fulfilled, one just has to take the gravity mediated effects, Eq. (|2]), into account 
and the behaviour becomes of the gravity mediated type at the time of condensate 
formation. 

We could also include thermal corrections to the potential, Eq. ([!]), as has been 
done recently J7J [T2|, [H5|, £2IJ. Their effect is of the same form as Eqs. (|2],|3]) if one 



replaces either m 3 /2 or m$ by the temperature T ^> TO3/2, rn<$>. In principle one would 
have to include also the gravitational corrections in gauge mediated case, too, but 
since all the relevant quantities of the AD condensate are determined by the most 



dominant mass term we only have to pay attention to one mass term at a time in the 
simulations. 

The main contribution to the evolution of the flat direction comes from the homo- 
geneous mode whose equation of motion is 

dV 

$ + 3H$ + — - = 0. (4) 

It is easier to analyze the behaviour of the equations of motion by parameterizing 
the field by $ = (pe td , for numerical purposes in Section 3 the parameterization 
$ = -^(01 + i<fa) is used. Then one obtains 

dV 

<p + 3H(f) - 9 2 <p + — = 0, 



In order to understand qualitatively the behaviour of the equations of motion 
Eq. (JD we find the maxima and minima of the potential, which are achieved with the 
phase values 

a = _I ow*™ ( ' A ' m3 / 2 sin ^ + ^ ± W H M0a + Ox) \ 2n + b 

d aiClan {\A\m 3/2 cos(9 A + e x ) + \a\Hcos(9 a + 9 x )) + d ^ W 

where 6 = for a maximum and 6=1 for a minimum.^ This can be simplified for 
H ^> m 3 / 2 for which the phase minimum is at d9 = —(9 a + 9\) + (2n + 1)tt, while for 
H <C itl 3 /2 the minimum is at d9 = —(9a + 9\) + (2n + 1)tt. For H ^> m 3 / 2 , m$ one 
can solve the minimum for cf) approximately as 



V2 



ilf d " 3 H 



p 



2|A|(d-l) 



f(9) + Jf(9y + A(d-l)c H 



l/(d-2) 



(7) 



where 



f(9) = \A\^cos(d9 + 9 A + 9 x ) + \a\cos(d9 + 9 a + 9 x ). (8) 



H 

These will be needed for numerical simulations as the initial conditions. 



2 Note that this is valid only if the denumerator is positive. If the denumerator is negative one has 
to switch 6 = for a minimum and 6=1 for a maximum. 
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Now we can analyze the AD condensate formation qualitatively by noting that the 
AD potential is analogous to potentials giving rise to a first order phase transition 
in the direction and in the F-term case there is also a second order transition in 
the 9 direction. However, in this case there appears no tunneling. In this case the 
Hubble parameter determines the behaviour of the potential, not the temperature as 
usual. During inflation the AD field $ settles to one of the symmetry breaking minima 
in Eqs. (§,0) in the F-term case. In the D-term case the phase 9 stays random and 
only acquires a minimum value Eq. (|7]) [II]. After inflation the Universe becomes 



matter dominated by virtue of inflaton oscillations and the Hubble parameter reads 
H = 2/(3i) so that the minimum, Eq. fl7|), evolves towards a lower energy scale. <fi 
evolves close to a fixed point value of the equations of motion Eq. (|5|) |TI|] , which is 
slightly larger than the minimum (see Eq. (]?])), and given by 

Eventually the Hubble mass term becomes equal to the mass term coming from 
soft SUSY breaking, Eqs. (@,||), = V m {§)/\§\ 2 . At this scale the phase transition 
from <p = m i n , Eq. (|7|), to <p = occurs and the AD field $ starts to rotate in the 
pit of the symmetry breaking minimum. For the gravity mediated case the phase 
transition happens when H pt ~ m 3 /2- In the gauge mediated case it happens at 
Hpt ~ 2. . .5m§\\\ 1 /( d ~ 1 \m$/M p p d ~ 3 ')/( d ~ 1 \ which can be calculated by noting that 
the logarithm is on the range 10 . . . 100 for a large range of and m$. After a while 
the symmetry breaking minimum vanishes while the AD field continues its rotation in 
a spiraling orbit around $ = 0. This is what charges up the condensate, the charge 
density of which is defined by 

q = - $$*) = 9<p 2 . (10) 

i 

The charge density in the co- moving volume asymptotes to a constant when H ~ 
O.lf/pf. Fig. |] shows the condensate formation. More on this in Section 3. 

The subsequent behaviour is determined by the mass term of the field In the 
original AD scenario [|T], the mass term was given by V m (<&) = m 2 |$| 2 and the 



condensate decayed into $-quanta. The evolution is however totally different if the 
potential grows slower than $ 2 as is the case considered in this paper. Then the 
condensate fragments into non-topoligical colitons called Q-balls || ||. 



2.2 Analytical properties of the AD condensate 

The charge density, Eq. ([T0|),is calculated using the second equation of Eq. (|D in the 
form d/dt(90 2 R 3 ) = -R 3 dV/d6, where R = R (t/t ) 2 / 3 is the scale factor of the 
Universe, giving approximately 

dV A ( M d p ~ 3 H \^ 
g ~W~ m3/2 { |Aj J ' ( ' 

where H ~ H pt at the time of condensate formation. 

When the condensate has been charged, we can approximately solve the equations 
of motion by neglecting all the non-renormalizable terms and the Hubble mass terms 
in Eq. ([]]) and retain only the mass term in the equations of motion Eq. ([5]): 

4> + 3H<p - e 2 (f) + c^' 1 = o 

0+ \3H + ) 6 = 0. (12) 



The second equation of Eq. dl2[) merely gives conservation of charge density in the 
co-moving volume Q = 6<f) 2 R 3 = const. In the first equation of Eq. ( |T2"D the derivative 
of the mass term has been approximated by C<p n ~ l (for the potential V(<p) ~ <f) n ). 
For the gravity mediated case C = 2^ K ^M 2 ^ K ^rriy 2 and n = 2 — 2\K\; for the gauge 
mediated case C = 2m% and n = can be used since the derivative of the potential 
has this behaviour for large 0. There then exists a fixed point solution (neglecting 
terms of order H 2 which is consistent since these terms there neglected in Eq. (|T2"D, 
too) with 

0= (l")" +2 R '^- (13) 

This behaviour is exactly the same as was obtained for coherently oscillating conden- 
sates in |16| . Eq. (p~3|) gives 4> ~ R~ 3 /( 2 ~\ K \) for the gravity mediated case and 4> ~ R~ 3 
for the gauge mediated case. The usual method of analyzing the stability of the fixed 



point by calculating the eigenvalues of the matrix of the linearized equations does not 
work, because Eq. ([12]) are non-autonomous differential equations. 

Another feature of coherent condensates with potential growing slower than $ 2 is 



that the condensate has a negative pressure f|, [16], |2j| . As a consequence it becomes 
unstable since the energy in a co-moving volume increases, as can be seen from the 
continuity equation 

dV 

p + 3H(p + p)- — = 0, (14) 

where p = |$| 2 — V(<&) and p = |$| 2 + V(§) and the partial derivative of potential with 
respect to time appears, because the potential depends explicitly on time through H. 
However, after condensate formation that term is negligible, since H 2 <C K($)/|<1 ) | 2 . 
The equation of state for coherent condensates can be found by averaging w = p/ p = 
2|$| 2 /p — 1 over a cycle of rotation as done for oscillating condensates in [TO], so we 
want to calculate the time average of J £ dt |$| 2 /p @. The calculation can be done 
by noting that \<&\dt = c/|$| is the arc length in the $ space. Now we get for w 

w = 2° — tt?-!, (15) 



c 



where we are now integrating over the orbit of the AD-field <£>. This is exactly the 
same as was the case for pure oscillation in [OJ except that now the integration is over 



a different orbit and the energy density is equal to p = V max + V m i n (for pure oscillation 
Knin = 0). Usually the orbit is quite close to an elliptical orbit with a semimajor axis 
A and a semiminor axis B giving the arc length as 



d|$| 



B 2 ip 2 dip i . , 

1 (16) 



and the integration limits are < ip\ < A (B < A). Eq. (15) cannot be given in a 
closed form except in a few cases: for circular orbit B — A from which follows that 
Vmax = Vmin and p = 2V , giving the result w = 0; for pure oscillation B = (no charge 



in the condensate in this case) when V($) ~ |$| n , giving w — (n — 2)/(n + 2) [16 
and if n = 2 for all values of < B < A, we get w — 0. 
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Figure 1: Pressure-to-energy ratio, w, plots (a) gravity mediated case for K = —0.1 (solid) 
and K = —0.01 (dashed), (b) gauge mediated case plotted against D = <t>max/ (^ m %) f° r 
e = 0, 10 -8 , 10 -4 , 0.1 with solid, dashed, dash-dotted and dotted lines. 



Let us generalize the result of |L6[ by considering elliptical orbits in a potential 
V($) ~ |$| n . We can simplify Eq. fllS]) to depend only on the ratio of semi-major 
axes e = B/A (actually e = yl — e 2 where e is the eccentricity of the ellipse) 



fdx{l + i 



W 



2 (l-x 2 ))' 



/2H1/2 



1/2 



1 - 



(x 2 +e 2 (l-x 2 )) n/2 



-1/2 



(17) 



where a change of variable tpi = Ax was done. When e goes from to 1, the orbit 
changes from pure oscillation to circular rotation. At the same w changes from (n — 
2)/(n + 2) to zero. Thus for n < 2 the pressure is negative. We have plotted in 
Fig. 0(a) the result of the numerical integration of Eq. ( |17D in the gravity mediated 



case. For pure oscillation the pressure has been calculated to be p ~ yp PHI , which is 



negative when K < 0. 

For the logarithmic mass term of the gauge mediated scenario, Eq. ([3D, the calcu- 
lation of pressure is more involved. If we approximate the potential as a constant, as 
was done in ]3[], we get for pure oscillation w = — 1 by putting n = 0, thus yielding 
the smallest value of w possible. With the logarithmic mass term we can evaluate w 
numerically by making the approximation that the orbit is elliptical. (We shall see in 
Section 3 that the orbit of the AD field $ is not really elliptical but looks more like a 
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rotating trefoil.) Thus we obtain 

r .2 \ 1/2 



lag(l+r>(x 2 +e 2 (l-a! 2 )))] 1/ ' 2 
log(l+-D)+log(l+De 2 ) 



w = 2^ 1, 

/ „ n s 1/2 r ,/,,„/.,, i,-, oxwl— 1/2 



/ f 1 + f 

o V x ~ 



log(l+D{x 2 +e 2 (l~x 2 ))) 
log(l+D)+log(l+De 2 ) 



-1/2 



where D = A 2 /(2m|). In Fig. [I] we plot the numerical results of Eq. (|T8| ) for several 
values -D. 

We have seen that for potentials that grow slower than |$ 2 | , the pressure is negative 
unless the orbit of the AD field is circular. In an expanding Universe a circular orbit for 
a spatially homogeneous condensate is impossible because of the dissipation. Therefore 
the energy density of the AD condensate in the co-moving volume increases. Thus the 
condensate is unstable and the most unstable orbit is purely oscillating orbit, where 
the charge of the condensate is zero. This feature was also seen in the simulation of 
Kasuya and Kawasaki ||, where they showed that there exist more unstable modes 
for pure oscillation compared to the circular orbit. This is the instability that causes 
the condensate to fragment. While it would be natural to expect that the condensate 
decays into $-quanta, for potentials growing slower than $ 2 the minimum energy 
configuration is an ensemble of Q-balls. 

2.3 Q-balls 

Q-balls have the generic form 

$(x,t) = i=0(r)e^ (19) 

where <f>(r) is a decreasing function of r. The value of rotation velocity u and the form 
of 0(r) depend on details of the potential. Note that a Q-ball is rotating circularly 
in the field space. Unlike for the coherent condensate, this is now possible because 
the Q-ball is spatially inhomogeneous. The different kinds of Q-balls considered in 
the literature can either have a narrow well-defined edge, in which case they are called 
thin- wall Q-balls or their boundaries are not localized, in which case they are called 



thick-wall Q-balls [22]. In this paper we concentrate only on thick-wall Q-balls. 
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Thick-wall Q-balls arise in supersymmetric theories with gauge and gravity me- 

[. In the gauge mediated case the field value 



diated supersymmetry breaking |23 



inside the Q-ball, energy and radius are are given by |p3 



m* i 



E 



Rc 



1 



\/2m, 



(20) 



where the profile of the field is well approximated by a kink-solution po[ . In the 
gravity mediated case the corresponding quantities are |E3 





r 2 


>(r) 


= 4>oe ^ 


E 


= rn 3/2 Q 




1 


Rq 




J\K\m 



(21) 



3/2 



There is also the possibility having mixed Q-balls with both mass terms Eqs. (^ |3|) of 
gravity and gauge mediated cases together p6| . 



The energy (mass) of the Q-ball has to be less than the mass of Q free scalars 
i.e. E < m\Q\ in order to be stable with respect to decay into free scalars. For this 
reason the AD condensate, which is a spatially homogeneous rotating field, fragments 
into Q-balls and not to free scalars. The interesting cases for us are the thick-walled 
Q-balls. 



2.4 Condition for thermal equilibrium 

In || it was noticed that for large values of energy-to-charge ratio x (or, to be more 
precise, for large \x\) the end product of condensate fragmentation is an ensemble of 
both Q-balls and anti-Q-balls. The total charge they carry is much larger than the 
total charge of the condensate: Q + + Q_ 3> Qtot, where Q + (Q-) is the total charge 
carried by Q-balls (anti-Q-balls). This can be understood easily if both the total 
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energy and charge are stored in Q-balls, as has been argued in || [7], || |9], |10|] : 

£ to4 = J dpdQ(N+(p, Q) + AL(p, Q))^Jp 2 + m 2 Q = E+ + E- 

Q tot = JdpdQ(N + (p,Q) + N_(p,Q))Q = Q + -Q_, (22) 

where mq is given by Eqs. (|20|,pl|) and iV±(p, Q) are the number distribution functions 
of Q-balls, which at this point are arbitrary. We approximate the integrals in Eq. (E 



by replacing Q with an average charge Q > of a Q-ball and yp 2 + mq with an 
average energy E to obtain pa E(N + + AL) and Q± ~ Q{N + — iV_), where 
-W± = / dpdQN±(p,Q). Then we can approximate the energy-to-charge ratio x by 

mg - Q- 



The factor E / wig can be interpreted as 7 = 1 / yl — f 2 . For the gravity mediated 
case, Eq. (EH]), we see that if x 3> 1 then either 7 ^> 1, which would make the Q-balls 
ultra-relativistic, or Q + + Q_ 3> — Q_ = Qtot, which indicates that the total 
charges of Q-balls and anti-Q-balls are much larger than the original net condensate 
charge. Since it is unlikely that massive particles with m ~ 100 GeV . . . 100 TeV and 
Q ^> 1 would be ultra-relativistic, there has to be a large charge asymmetry in Q-balls 
for x ^> 1. For the gauge mediated case Eq. (|20|) a large number of Q-balls and anti- 
Q-balls can be produced even if x < 1, because Eq. ([23|) contains a numerical factor 
rriQ/mQ < 1 for Q > 1. 

Since the number of Q-balls and anti-Q-balls is large, it is likely that there are 
numerous Q-ball collisions, which rapidly thermalize the distibution. The condition 
for this to happen is that the collision rate of Q-balls is larger than the Hubble rate, 

T = n tot av >H = ^ (24) 

where n to t is the number density of Q-balls, a ~ tcRq is the cross-section of a Q-ball 
collision and v is the average velocity of a Q-ball. Using Eq. (p3|) we can write the 
rate as 

rPax — q tot 7rR 2 n-, (25) 
mn v 7 
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where q to t = Qo{to/t) 2 is the total charge density with q the initial charge density 
of the condensate at time to- in order to get the most conservative bound for the 
collision rate we approximate the average charge of the Q-ball, Q, with the maximum 
charge of a Q-ball, Q max - Fragmentation of a d — 4 condensate was considered by 
Kasuya and Kawasaki both in the gravity mediated || and gauge mediated case 
H They performed numerical simulations which begun at to = 2/(3m 3 / 2 ) (grav.) 
and at to = \/20o/(3m|) (gauge) with initial charge density go = m 3/20o anc ^ f° un d 
that the maximum charge of the Q-balls formed is Q max ~ 6 • 10^ 3 4>lm^ 2 (grav.) and 
Qmax ~ 6 ■ lCT 4 0QmJ 4 (gauge). The Q-balls are formed at t/ ~ 5 • 10 3 m^ 2 (grav.) 
and tf ~ 5 ■ lO 5 ^,^ 1 (gauge) |§. Now Eq. fl24]) becomes 



j 1 Grav. med. 

x-> I / \ (26) 

7 io-2(Mev) Gauge med. 



™3/2 

where we adopted \K\ = 0.1 and O = (M^~' 3 H y|A|) 1/(d ~ 2) . Eq. (H) shows that for 
large values of x the condition of thermal equilibrium is fulfilled without assuming any 
specific form of the distribution. If the Q-balls have a relativistic average velocity, 
v ~ 0.1 ... 1, then even with x ~ 1 (grav.) and x < 1 (gauge) Eq. (^) is fulfilled. 

One should note that actually the condensate is not formed strictly at the initial 
times used in Eq. (f2T)|), which were taken from || [7|, §], but a little later, as can been 
seen from simulations in Section 3. For the gravity mediated case the formation time 
is approximately t ~ 10. . . lOOm^ and for the gauge mediated case t ~ lO 5 /?^ 1 . 
Thus the right-hand side of Eq. ( PB"D is even smaller by a factor of 0.01. ..0.1. 

Another matter is the value of the average charge Q. If we use a thermal distribu- 
tion we get from |J that Q < 10 1 in the gravity mediated case. This is orders of mag- 
nitude smaller than Q max so that the condition of thermal equilibrium, Eqs. (p^, p6]), 
is fulfilled self-consistently. It was also shown in [[| that even with a low energy-to- 
charge ratio such as x ~ 10 the thermal distribution is a very good approximation to 
the true distribution. In the gauge mediated case therefore the issue of thermalization 
of the distribution has not yet been checked in the existing simulations. The existence 
of thermal distributions has been established in the gravity mediated case for both 2 
and 3 spatial dimensions |5|, [TIJ . 
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The average velocity, v, is undetermined from these arguments. Energetically the 
existence of relativistic Q-balls is possible, since even with v ~ 0.1 we obtain 7 ~ 1.005, 



leaving Eq. (|2~3"D essentially unaffected. Only if v ~ 1 we obtain 7> 1. Therefore it is 
natural to expect that Q-balls are relativistic. 

We can calculate the absolute minimum of \x\ by treating 0, <fr, 9 and 9 as inde- 
pendent variables to obtain 



■ 1 p \4> 2 + \o 2 4> 2 + v{^e) > / 2V(0,g min ) 

m\q\ ~ m\9W ~ V ™¥ 2 ' 



where we can approximate the potential, V, by its mass term alone. In the gravity 
mediated case this gives |x| > 1.01 (for 0/M = 0.1 and K = —0.01), and in the 
gauge mediated I > 10 3 (10 8 ) for d = 4 (6). In practice these are not realized. 

Because minimum |x| would require a circular orbit. In the next Section we will discuss 
the realistic values of x. 



3 Numerical simulations 

3.1 Details of the numerical simulations 

We have simulated the evolution of the flat direction starting from the scale of infla- 
tion to the time of AD condensate formation. The numerical simulations were done 
by using routines for ordinary differential equations from the NAG library. We pa- 
rameterized the AD field $ as real and imaginary parts, rescaled the fields and time 
to be dimensionless and used logarithmic time in order to handle the scale difference 
between the scale of inflation and AD condensate formation 

1 / M d - 3 H\ 1/{d ~ 2) 

$ = 7^(01 + #2) , <]>i = f " |A| j X, , z = log(mt), (28) 

where m = 717.3/2 in the gravity mediated case and m = mq> in the gauge mediated 
case. With a redefinition of the phase of the field $ we put 9a = and 9\ = it so that 
the minima are at b = in Eq. (§) i.e. d# min = —9 a for large H, when a ~ 1, and at 
c^min = for small H. The initial time is the end of inflation so Zi = \og(2m/3Hj), 
where Hi = 10 12 GeV was chosen. Because of the fixed point behaviour Eq. (|9]), the 
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exact value of initial time is not relevant as long as the time, when rotation starts, is 
much larger than the initial time t rot ^> U. Then the equations of motion, Eq. 
become 



xi + Ci*i - c 2X i + Kxi + c 3 (xl + xl) {d - 2)/2 xi 
c,(xi + xi) {d - 1)/2 



A|m3/2 e z cos((d - 1)0) + ^! cos ((d _ i)g + e a ) 



ni 



X2 + C1X2 - C2X2 + V^X2 + C 3 (xi + X 2 



2\(d-2)/2 



X2 + 



^4(X? + ^) (d - 1)/2 

where 



|A|m3/2 sin((d - 1)6) + sin((d - 1)6 + O ) 



??? 



0,(29) 



c 2 
c 3 

^4 



d-4 
d-2 

d-3 
(d-2) 2 + 

d-1 

9 ■ 2 d - 4 
1 

3 . 2(<i-4)/2 



4c H 



(30) 



and 



.-,2- 



1 - |K| log 



V 



,2z 



1 (M^m^2 

1 H — < 

2ml V IAI 3 



' M^m 3/2 2\ 2/{d - 2) X l + X t 
|A| 3/ 2M 2 

(xl + xl) 



3 / 

2/(d-2) 



(Grav.med.), 



(Gauge med.) (31) 



and = arctan(%2/xi)) c # = 1^1 = 1? I a l = 0, 1. We chose the renormalization scale 

/ ._„ \l/(d-2) 

in the gravity mediated case to be M = [M^m 3/ 2/\M) as was done in ||18||. 

However, we have also checked the results for larger values of M that have been used in 
simulations all the way up to Planck scale |7|, §, ||, [l(| . We comment on the differences 
when discussing the results. The initial conditions Eqs. (0,0) in these variables read 

l/(d-2) 



X1+X2 



V2 



a\ + sJ\a\ 2 + A(d- l)c H 
2(d- I) 



G [-*§], a~0 (D-term) 

a ~ 1, 6 a e [-7T, tt] (F-term) 



(32) 
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Figure 2: Affleck-Dine condensate formation (a) gravity mediated case with d = 4 (solid) 
and d = 6 (dashed) and (b) gauge mediated case with d = 4, m$ = 1 TeV (solid) and 
m<j> = 10 TeV (dashed) with initial condition 6% = — 7r/10 and a = in both figures. 

We have restricted the initial phase, 9i to this sector since changing the initial phase 
by 6 — > 6 + ^jH gives the same result for charge and energy of the condensate. The 
solution of $ in another sector of initial conditions is achieved by a phase shift of 
2ixn/d from the sector prescribed. 

3.2 Outline of the numerical results 

Fig. |] shows examples of the rotation of the AD condensate. In the gravity mediated 
case, Fig. 0(a), we see that the orbit is a spiraling ellipse and in the gauge mediated 
case, Fig. |2|(b), a precessing trefoil. From Fig. ^|one can see that there is a twist on the 
orbit much before the rotation starts properly. This is the time of the phase transition, 
when the AD field $ starts to rotate in the pit of the symmetry breaking minimum. 
The rotation begins when the symmetry breaking minimum is the vacuum, ends when 
it has become a false vacuum and twists when the false vacuum has completely vanished 
forming a kink on the orbit. It is possible to produce a condensate through a second 
order phase transition, too, but the charge in that case would be small. It should 
also be noted that in the gravity mediated case condensate formation starts when 
chH 2 ~ m^ 2 for all values of d, A and a. In the gauge mediated case the condensate 
formation starts at ChH 2 ~ m|/|$| 2 , so that the formation happens earlier if the mass, 
m$, is increased, as can be seen from the different positions of the kink in Fig. 0(b). 
The quantities of interest are the charge density q, the energy-to-charge ratio x, 
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the pressure-to-energy ratio w, the ellipticity of the orbit e and the decay of the 
AD field amplitude <fi. The charge density contributes to the baryon-to-entropy ratio 
and is depicted in Figs. |3] and |7|. The energy-to-charge ratio x, which gives the 
relative amount of Q-balls, anti-Q-balls and their kinetic energies, Eq. (|23|), is shown 
in Figs. |] and |8]. The pressure-to-energy ratio w, which gives the equation of state 
and indicates whether the condensate is unstable or not, is shown in Figs. § and |9|, 
where we also show the ellipticity of the orbit e (the calculation of the average w is 
technically somewhat problematic, because of the rapid oscillations). The of the AD 
field amplitude 0, to verify the proposed behaviour Eq. (0), is shown in Fig. || 

We display the results for both D- and F-term inflation in the gravity and gauge 
mediated cases: we adopt dimensions d = 4, 5, 6, 7 (gravity) and d = 4, 6 (gauge). 
The gravitino mass in the gauge mediated case is chosen to be 777,3/2 = 10~ 5 m<i, (d = 4) 
and 7773/2 = 10 _9 777$ (d — 6). The variation in 7773/2 only results in variation in the 
charge density and therefore affects x inversely. In the gravity mediated case the 
oscillation of x asymptotes at around t = lOOmT^. In the gauge mediated case there 
is no asymptotic behaviour but x continues to increase steadily. Therefore in the gauge 
mediated case the quantities in the Figures are presented for different times. 

We have summarized the values of q, x, w and e in Table 1 (gravity) and Table 2 
(gauge). 

3.2.1 Gravity mediated case, D-term 

In this case a = and the mass term is given by Eq. (0). Initially the phase, 8i, 
is random and is at the instantaneous minimum Eq. ([?]). The parameter values 
adopted are M = (M p d - 3 m 3/2 /|A|) 1/(d - 2) ; \A\ = c H = 1; K = -0.01 (and K = -0.1); 
and d = 4, 5, 6, 7. 

In Fig. |3|(a) we show the time evolution of the co-moving charge density, q(R/ Ro) 3 - 
One can see that it becomes a constant when H ~ 0.1 . . . O.OI7773/2. In Fig. |3](b) 
the charge density in the co-moving volume is plotted against the different initial 
conditions. It can be seen that the charge density is positive for — 7r < ddi < 
and negative for < d9i < ir. However, this behaviour depends on the values of 
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Figure 3: Charge density in the co- moving volume in the gravity mediated case vs. 
time in logarithmic units of time, (b) D-term (a = 0) and (c) F-term (a = 1) case 



(a) 
with 



d = 4, 5, 6, 7 (solid, dash-dot, dashed and dotted lines), K = —0.01 and t = 100m 



3/2- 



ch and A. For instance, here we have chosen ch = 1 whereas in 11 the choice 



c h — 9/4 was made resulting in a negative charge density for some values d9i < 0. 
This comes about because the AD-field $ manages to rotate around the symmetry 
breaking minimum before it vanishes, thereby gaining phase motion in the opposite 
direction (we have verified this behaviour). Therefore the sign of the charge cannot be 
given generically for any range of initial conditions. Only the magnitude of the charge 
density can be given and is about the order depicted in Fig. 0(b). Typically the charge 
density lies in the range q(R/R ) 3 ~ 0.1 . . . 0.6m 3 / 2 {Mp~ 3 m 3 / 2 /\\\) 2 ^ d ~ 2 ^ ■ This results 
in a charge-to-entropy ratio q/s ~ 0.01 ... 0.1 (T R /m s / 2 ) (m 3 /2/|A|M p ) 2// ^~ 2 \ where 
Tr is the reheating temperature. Here we have assumed that the entropy density 
s « 4pi/Tr, where the inflaton energy density is pi ~ 3M 2 H 2 . For d — 4 and 
= 100 GeV this would require Tr ~ 10 9 . . . 10 10 GeV if the baryon-to-entropy 
ratio, tib/ s ~ 10~ 10 , is to be explained through AD mechanism. For d = 6 we would 



obtain Tr ~ 10 . . . 100 GeV. These are of the same order as the results of JTT[] such 
that the radiative correction does not produce any significant variation. 

In Fig. f|(a) we depict the time development of the energy-to-charge ratio x. One 
can see that x oscillates strongly with a decaying amplitude and the overall behaviour 
is that |x| increases. At t — lOOm^ the oscillation amplitude has dampened enough 
such that x is approximately a constant. In Fig. |](b) we plot x against the different 
initial conditions. There one can see that as d increases, \x\ also increases. Roughly 
one third of the total range of the initial conditions were found to have \x\ > 10 
while two thirds have 1.1 < \x\ < 10 (lower limit is the minimum value for d = 4). 
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Figure 4: Energy-to-charge ratio, x, in the gravity mediated case vs. (a) time in logarithmic 
units for d = 4, 6, (b) D-term (a = 0) and (c) F-term (a = 1) case with d = 4, 5, 6, 7 (solid, 
dash-dot, dashed and dotted lines), K = —0.01 and t = lOOm^. 

Therefore more simulations for \x\ < 10 should be made to establish whether the 
thermal distribution is valid or not. If \K\ or M are increased from the values used 
here, \x\ becomes slightly larger. 

The oscillation of pressure-to-energy density ratio, w, is plotted in Fig. ||(a), which 
shows that the oscillation of x corresponds to oscillation of pressure from positive to 
negative. The average pressure is slightly to the negative side, resulting in an increasing 
x. We have calculated the average w in Fig. |5|(b) at t ~ lOOm^ with different initial 
conditions. In Fig. |5](c) the ellipticity of the orbit, e, is plotted to show that w is 
more negative if e is small. However, it should be noted that w achieves values which 
are more negative than the absolute lower bound coming from pure oscillation, see 
discussion after Eq. ([[]]). Therefore we also show the average w at t — 300m^ 2 in the 
same figure. This behaviour might be due to dV/dt in Eq. (0), whose contribution 
decreases as time evolves. It is also due to the rapid oscillation of w that there are 
numerical inaccuracies in our calculations. However, we believe that this does not 
affect our main conclusions. 

From Fig. |](a) one can see that the amplitude, 0, decays as (R/ Rq)~ 3 ^ 2+k ^ where 
K = —0.01. We checked that for K = —0.1 the decay of <p tends to slow down. 
This is due to the fact that for K = —0.1 the approximation of the mass term as a 
polynomial breaks down and in effect gives rise to smaller \K\. However, the effect is 
not significant unless M very large. 
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Figure 5: Pressure-to-energy density ratio, w, in the gravity mediated case vs. (a) time in 
logarithmic units for d = 4 and (b) different initial conditions for d = 4, 6 and (c) ellipticity, 
e, vs. initial conditions, d = 4 (thin lines), d = 6 (thick lines), D-term (solid), F-term 
(dashed) with K = —0.01 and t = 100m~^ 2 . In plot (b) there is plotted w at t = 300m^^ 2 
with dotted lines for d = 4 D-term case. 

3.2.2 Gravity mediated case, F-term 

In this case \a\ — 1 and the initial conditions are given by Eq. (^). There is however no 
practical difference with respect to the D-term case, only minor quantitative changes. 
In Figs. (|3],^|)(c) we have plotted the charge density in the co- moving volume and x 
to compare them with the D-term case. Charge densities are approximately similar 
for d = 5, 6, 7 but for d = 4 there is a difference. This is probably due to the nature 
of the fixed point Eq. (|SJ): for d = 4 it is marginal whereas for d > 4 it is attractive 



llj. Now \x\ is slightly larger than in the D-term case, the minimum being \x\ > 2. 
The orbits are also slightly more elliptic, as can be seen in Fig. |5|(c) resulting in more 
negative pressure, see Fig. |](b). 

Table 1. Gravity mediated case 



d 


4 


6 


q(R/Ro) 3 / m 3/2 (M«- 3 m 3/2 / \ A | f'^ 


~ 0.1... 0.6 


-0.1. ..0.3 


\x\ 


> 1.1 


> 2.5 


w 


(0.05 . . . 1)1^1/2 


~ -(0.2... 1)1^1/2 


e 


< 0.5 


< 0.2 



19 




Figure 6: Time development of the field amplitude (a) 4>(R/Rq) 3 ^ 2+k ^ in the gravity me- 
diated case for d = 4 and K = —0.01; gauge mediated case 4>(R/Ro) 3 for (b) d = 4 and (c) 
d = 6. 

3.2.3 Gauge mediated case, D-term 

In the gauge mediated case the mass term is given by Eq. (||). Now a = 0. We 
have checked the cases d = 4, 6 with 7713/2 = 10~ 5 , 10 _9 m<j>, \A\ = ch = 1 and 
m$ = 1, 10, 100 TeV. Here we have chosen 7113/2 such that a maximal amount of 
charge, and therefore the minimum value of \x\ is produced. If 1713/2 were larger, we 
would have to include the mass term of the gravity mediated case and the behaviour 
would effectively be that of the gravity mediated type already discussed. 

In Fig. |^(a) the time development of the charge density in the co-moving volume, 
q(R/Ro) 3 , is shown for d = 4. One can see that it reaches its asymptote earlier for 
large m$. The charge in dimensionless units is different for different m$. This is due to 
the fact that one cannot simultaneously remove the dependence on m$ and 7713/2 ■ The 
resulting charge density in the co-moving volume is q(R/Ro) 3 ~ 0.1M p m|/|A| for d = 4 
and q(R/ Ro) 3 ~ 1000m$(Mpm$/|A|) 1 / 2 for d = 6. This results into charge-to-entropy 
ratios q/s ~ 0.1T R /M p for d = 4 and q/s ~ IOOTrM" 1 / 2 ™^ 172 for d = 6. 

In Fig. |8](a) the time development of energy-to-charge ratio, x, is shown for d = 4. 
One can see that there is no oscillation in the gauge mediated case in contrast to the 
gravity mediated case, x increases faster, which is due to the behaviour <p ~ -R -3 , see 
Fig. [|(b), which causes the rotation velocity 9 ~ R 3 to increase rapidly. In Figs. |8](b) 
and (c) we show the d = 4 and d = 6 cases. One can see that for d = 4 one obtains 
\x\ > 10~ 2 whereas d = 6 results in \x\ > 10~ 5 (for m$ = 1 TeV, which gives the 
lowest bounds). If 777,3/2 decreases, the charge density becomes smaller, thus making 
which is inversely proportional to charge density, larger. In a realistic case one 
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Figure 7: Charge density in the co-moving volume in the gauge mediated case vs. (a) time 
in logarithmic units of time d = 4, (b) d = 4 and (c) d = 6 with D-term (thin lines, a = 0) 
and F-term (thick lines, a = 1) and m$ = 1, 10, 100 TeV (solid, dashed, dotted lines) at 
t = 4 • 10 5 mj\ 10 5 m5\ 4 • lO 4 ™^ 1 (d = 4) and t = 4 • 10 9 , 10 9 , 4 • lO 8 ^ 1 (d = 6). 




lo 9,oCV> 



Figure 8: Energy-to-charge ratio, x, in the gauge mediated case vs. (a) time in logarithmic 
units of time d = 4, (b) d = 4 and (c) rf = 6 with D-term (thin lines, a = 0) and F- 
term (thick lines, a = 1) and m$ = 1, 10, 100 TeV (solid, dashed, dotted lines) at t = 
4 • 10 5 m5\ lO^J 1 , 4 • lO 4 ^ 1 (d = 4) and t = 4 • 10 9 , 10 9 , 4 • lO 8 ™^ 1 (d = 6). 

should take the gravity mediated mass term, Eq. (|]), together with the gauge mediated 
mass term, Eq. ([3|), for these borderline cases. This would increase \x\, so that the 
true minimum values are likely to be an order of magnitude larger. One should also 
note that in order for the AD condensate to form via the gauge mediation mechanism 
one requires m 3 / 2 < 0.1 ... 10 GeV for d = 4 or m 3 / 2 < 1 . . . 100 eV for d = 6. For 
d = 4 condensate formation is quite possible, but for d = 6 it is likely that a gravity 
mediated mass term must be included. Then it would appear that the AD condensate 
forms through gravity mediation, but the Q-balls via gauge mediated mechanism or 
via a mixture of these |26| . 

In Figs. |](a) and (b) we show the time development of the pressure-to-energy 
density ratio, w, for d = 4 and d = 6. One can see that the pressure is always negative. 
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Figure 9: Pressure-to-energy density ratio, w, in the gauge mediated D-term case (o = 0) 
vs. time in logarithmic units for (a) d = 4 and (b) d = 6. (c) plot of ellipicity of the 
orbit where d = 4 (thin lines) and d = 6 (thick lines). The values of the scalar mass where 
m<j> = 1, 10, 100 TeV with solid, dotted and dashed lines. 



The calculation of average pressure is even more involved than in the gravity mediated 
case, since the oscillation frequency becomes very large. In Fig. P(c) we show the 
ellipticity as a function of different initial conditions. We note that quite generically 
e< 0.1. 

In Figs. P(b) and (c) we display the time development of the field amplitude, 
<f)(R/Ro) 3 , for d = 4 and d = 6. Note that <fi is oscillating, thus confirming the R~ 3 
behaviour derived analytically in Eq. (p~3|). 

3.2.4 Gauge mediated case, F-term 

Now \a\ = 1. This case is essentially the same as the D-term case, as can be gathered 
from Figs. [j], |8] and |9|. Note again that \x\ is slightly larger compared to the D-term 
case while the orbits are more elliptical, as in the gravity mediated case. 



Table 2. Gauge mediated case 



d 


4 


6 


q{R/R ) 3 / m $ (M p d - 3 m $ /|A|) 2 /( d - 2 ) 


~ 0.1 . . .0.6 


~ 1000 . . . 6000 


\x\ 


> io- 2 


> io- 5 


w 


0.9. . . - 0.8 


0.95... -0.9 


e 


< 0.15 


< 0.05 
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4 Conclusions 



We have studied the evolution of the flat direction field starting from the end of infla- 
tion up to the formation of a coherently rotating Affleck-Dine condensate. The original 
studies |], |TTJ used a mass term V m {§) = m 2 |$| 2 . We have extended this to radia- 
tively corrected mass terms both in the gravity and gauge mediated case. Typically in 
both of these cases the mass term grows slower than |$| 2 . This is the reason why the 
pressure of the condensate is negative, thus making it unstable with respect to spatial 
perturbations. The other necessary requirement for a negative pressure was found to 
be the ellipticity of the rotation of the AD field: for circular orbits the pressure of a 
coherently rotating field is zero regardless of the potential. Negative pressure is largest 
for simply oscillating field. Its charge would be zero, resulting in an infinite energy-to- 
charge ratio, which would make the oscillating configuration the most unstable one. 
This can also be seen from the analysis of Kasuya and Kawasaki || for the growth of 
perturbation modes: for orbits close to the circular there is only one band of growing 
modes, but for an oscillating one there appear several growing modes. 

We studied numerically the time evolution of various quantities, such as charge 
density q, energy-to-charge ratio x, pressure-to-energy density ratio w, ellipticity of 
the orbit e and the evolution of the AD field amplitude <p. As an initial condition the 
homogeneous field was chosen to be at rest at the ground state of the potential. Only 
the initial phase of the field was varied over all the possible values for both D-term 
and F-term inflation. We found that there is no qualitative difference between the two 
cases. We checked both the gravity mediated case and the gauge mediated case with 
the dimension of the non-renormalizable operators d = 4, 5, 6, 7. (For d = 4, K = 



our results agree with IjlTJ.) In contrast to (TT| the energy-to-charge ratio, x, was found 
not to asymptote to a constant value. Instead \x\ follows an increasing curve due to 
negative pressure with a minimum \x\ > 1.1 in the gravity mediated case, \x\ > 1CT 2 
for d = 4 and \x\ > 10~ 5 for d — 6 in the gauge mediated case. We should mention 
that the calculation of the average of w is plagued with technical difficulties. However, 
the qualitative behaviour of w is not dependent on these. We believe we have managed 
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to check numerically that the ellipticity of the orbit really affects the pressure exactly 
the way indicated by the analytical considerations. In addition the field amplitude 
was shown to decay _R~ 6 /( n+2 ) for potentials V(<f>) ~ (fi n , where n = would correspond 
to the gauge mediated case, which is in exact analogy to the coherent oscillation ||T6| . 



In conclusion, we have shown that the AD condensate is unstable due to negative 
pressure, which is caused by two effects: the potential growing slower than the field 
squared and the ellipticity of the orbit. We have also shown that after the fragmen- 
tation of the AD condensate the thermalization of the resulting Q-ball distribution 
depends on the energy-to-charge ratio of the condensate. In the gravity mediated case 
x was found to be such that thermalization is likely to be a generic feature, as sug- 
gested by analytical considerations. The gauge mediated case is complicated by the 
fact that condensate formation can take place while the field is located in the part of 
the potential dominated by gravity mediated mass term, whereas Q-balls can be cre- 
ated by virtue of gauge mediation. However, as discussed Sect. 3.2.3, thermalization 
appears to be very likely also in this case. 
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